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The formation of quantized vortices in trapped, gaseous Bose-Einstein con- 
densates is considered. The thermodynamic stability of vortex states and the 
essential role of the surface excitations as a route for vortex penetration into 
the condensate are discussed. Special focus is on finite-temperature effects 
of the vortex nucleation process. It is concluded that the critical angular fre- 
quencies for exciting surface modes with the relevant multipolarities yield, 
also at finite temperatures, the appropriate thresholds for the nucleation of 
vortices in dilute Bose-Einstein condensates, in fair agreement with the re- 
cent experiments. 
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1. INTRODUCTION 

The existence of quantized vorticity is an essential part of superfluid 
plienomena. Quantized flow with 2it phase winding around a nodal line of 
a coherent many-body wavefunction is what is usually understood by the 
concept of a simple vortex line in the framework of quantum mechanics. 
These topological defects have previously been created and observed, for 
example, in quantum fluids such as liquid ^He and '^He, superconductors, 
and most recently, starting in 1999, in the Bose-Einstein condensates of 
dilute alkali-atom gases.B'i 

In type-II superconductors, the order parameter is suppressed along 
quantized vortex filaments when the external magnetic field exceeds the 
lower critical field, Hd- In higher magnetic fields, the superconducting sam- 
ple is pierced by an Abrikosov lattice of vortices having normal cores. As 
the field strength is further increased, the superconducting state vanishes at 
the upper critical field, Hc2, as the average distance between the vortices 
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becomes comparable with the size of the vortex core. 

In superfluid hehum, vortices are usually generated by mechanically ro- 
tating the fluid-containing vessel. Formally, the angular rotation frequency 
for the helium liquids is the analogue of the magnetic field in superconduc- 
tors. When superfluid helium is rotated above 0,ci, the fluid is penetrated 
by quantized vortices. The determination of the exact value of fici for the 
helium superfluids is complicated, for example, by effects such as surface 
roughness of the vessel, remanent vorticity, barrier-tunneling events, and 
the difficulties in the imaging of the vortex cores. For higher rotation ve- 
locities, an increasing number of vortices enter the system forming a regular 
vortex lattice, but the predicted value for the upper critical frequency Qc2 
is unattainably high and the regime of such extreme rotation velocities is 
beyond experimental reach. 

In dilute Bose-Einstein condensate gases ,Bq3bQ'EI vortices have been 
generated by a variety of different methods. The first vortex was created 
by directly imprinting the characteristic 27r phase circulation in the two- 
component condensate.^ A frequently employed method for creating vortices 
and vortex lattices — with largest of them containing more than a one hun- 
dred vortices — has been the stirring of the condensate by a time-dependent 
trap anisotropy.iJ2l0'0'0llillil Vortex lattices have also been created through 
intrinsic rotation, in which case the thermal boson gas is first rotated by a 
trap anisotrox!ij_and the gas is subsequently cooled below the condensation 
temperature. In those experiments it is the normal gas surrounding the 
superfluid that defines the rotating environment for the condensate, instead 
of the externally rotated trap potential. Since the gaseous condensates are 
confined by magnetic fields there is no associated surface roughness as in the 
case of superfluid helium, and thus substantial differences in the nucleation 
processes are expected to occur between these superfluid systems. 

Moreover, vortices have been produced through the snake instability of 
decaying solitonsEj and in the wake of the turbulent flow generated by a 
moving object.0'lll Dark solitons have also been seen to decay into vortex 
rings.li3 Most recently, multiply quantized circulation has been created topo- 
logically in dilute Bose-Einstein condensates by adiabatically switching the 
direction of the magnetic field along the vortex axis .1^3 

The organization of this paper is as follows: In Sec. 2 we outline the 
formulation of the mean-field theory relevant for the subsequent discussion. 
The physical mechanism underlying the vortex-formation process is ana- 
lyzed both qualitatively and quantitatively at zero and finite temperatures 
in Sec. 3. Section 4 summarizes our conclusions. 
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2. STATIONARY MEAN-FIELD THEORY FORMALISM 

Our starting point is the second-quantized grand-canonical Hamiltonian 
7Yf7 = Ti. — CI ■ L — fiN defined in the frame rotating at an angular velocity 
n by 



fl2 

-— + VtrM - n + ihfl • r X V 
2m 



^(r) dr 



Wf7 = y ^^(r) 

+ l y ^'■f(r)^'1'(r)l'(r)^(r)dr (1) 

where g describes the strength of the two-body interactions between the 
atoms and Vtr{v) denotes the external potential used for spatially confining 
the atoms. 

The bosonic field operator is decomposed in the usual fashion as ^'(r) = 
^(r) -|- V'(r), where (p{r) stands for the macroscopic wavefunction of the 
condensate. The mean-field scheme corresponds to approximating the cubic 
and quartic operator products according toEj 

ijj^ (r) if) {r)ip{r) ~ m{r)ip^ (r) + 2n{r)Tp{r), (2) 

where we have introduced the noncondensate density n(r) = ('0^(r)'0(r)) 
and the anomalous correlation function m(r) = ('0(r)'0(r)) describing the 
condensate-induced correlations between the noncondensate particles. 

Inserting the Bogoliubov decomposition for the field operator and the 
mean-field approximations in Eq. (||) , we obtain an eff'ective quadratic Hamil- 
tonian, Ti-eff- The linear terms in Ti-eff are eliminated by the stationary 
generalized Gross-Pitaevskii equation, 

[£ - 5|0(r)n </.(r) + gm{r)<p^r) = 0, (3) 

for the macroscopic wavefunction (/)(r) and the chemical potential ^. Above, 

C = -— V2 + Vtr{v) -fi + 2g [|<^(r)|2 + n(r)] + i^fi • r x V, (4) 
2m ^ ^ 

and the condensate wavefunction is normalized to the fixed number of par- 
ticles, A^, according to / [|(/)(r)|^ -|- n(r)] dr = N. 

The remaining terms in the effective Hamiltonian may be diagonalized 
by using the canonical Bogoliubov transformation 



= [^iq(r)??q - v*^{r)r]l , (5) 
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and demanding the quasiparticle wavefunction amplitudes Mq(r) and Vq{r) 
to obey the coupled Hartree-Fock-Bogoliubov eigenvalue equations 



£uq(r) - MVdir) = ^qUq(r) 
£X(r) - 7W*nq(r) = -^q^;q(r) 



(6a) 
(6b) 



where Ai = g[m{r) + (^^(r)] and Eg label the quasiparticle eigenenergies. 
Moreover, conservation of the bosonic commutation relations implies the 
normalization condition 



[u*iir)uj{r) - v*{r)vj{r)] dr = 6, 



(7) 



for the quasiparticle amplitudes. The positive-norm quasiparticle states with 
negative eigenenergies Eg are referred to as the 'anomalous modes'. 

Within the mean-field approximation, the quasiparticles are noninter- 
acting, the corresponding creation and annihilation operators obeying the 
following equilibrium relations; 



iVirij) = ivWj) = 



(8a) 
(8b) 
(8c) 



where rii is the Bose-Einstein distribution function. Using these relation- 
ships, the self-consistent mean fields may be expressed in the form; 



n r 



m(r) 



Y,[n^{\u^{v)\^ + \V,{V)\^)+\V,{V)\ 
i 

- [2niUi{v)v*{v) + ni(r)<(r)] . 



(9a) 
(9b) 



Equations (^) , (^) and (^) constitute a self-consistent set of equations whose 
solution for a partially Bose-Einstein condensed system of atoms require 
iterative methods. 

The expectation value of the effective Hamiltonian We// is found to be 



(We//) 



2gh{r) 



(r) dr 



+ 



+ 



2|m(r)p (/)^(r)m*(r) (r)m(r) 



(Vi^(r)£^(r))dr. 



dr 



(10) 
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Taking Eqs. (|5|) and (^) into account, the last line in Eq.([To|) is expanded in 
the quasiparticle basis as 

(V'^(r)£V^(r)) = ^ [niV*,{r)Cui{v) + {m + l)vi{r)Cvl{v)] . (11) 

i 

Multiplying Eq. (^a) and the complex conjugate of Eq. ^o) by Ui{r) and 
Vi{r), respectively, we find 

(^t(r)£v3(r)) = ^[ni{Mu*{r)vi{r) + \ui{r)\^Ei) 

i 

+ im + l){Mviir)uUr) - \vi{r)fEi)] . 

(12) 

Within the framework of the mean-field approximation, Eqs. (P), the above 
expression reduces into 

(V^t(r)£^(r)) = J2 [n^{\u^{r)\^ - \v,{r)f) - \v,{r)\^] - Mm*{r). (13) 

i 

The quasiparticle normalization condition, Eq. (0), thus results after rear- 
rangement of the terms into 



/:-g^\cP\^-2gh{r) 



0*(r) 

+ ^niE,-^Ei j \vi{r)\^ dr. 



+ 



b{r) dr 
dr 



(14) 



The second line of Eq. ( |l^ vanishes since, by definition, (Ti-eff) is a real 
function, thereby ImJ ?7i*(r)(/)^(r) dr = 0, and we find 



(We//) 



r)| dr — 2g / \4>{r)\ h(r)dr — g / rh[r)(f)* (r) dr 



+ ^mEi-Y^Ei j \vi{r)\^dr 



(15) 



where the first line is expressed using the Gross-Pitaevskii equation. The 
entropy of the noncondensate atoms equals 



5" = -/cB^Klnnj - (1 -Frij) ln(l +ni)] 

i 

while the free energy of the gas is finally obtained from 

F = {neff)+lJiN -TS. 



(16) 



(17) 
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In a typical vortex-creation experiment, the Bose-Einstein condensate is 
rotated at a certain angular velocity by the time-dependent trap anisotropy. 
After the stirring, the condensate is allowed to equilibrate. Subsequently, the 
trap is turned off and the cloud is allowed to expand such that the possible 
presence of the vortex cores can be imaged. Computationally, vortices have 
been verified to enter the condensate from its boundary region by a number 
of simulations 

Characteristic for these rotating-bucket experiments has been the obser- 
vation that the vortices are only seen at substantially higher angular frequen- 
cies than what is expected from the basis of purely thermodynamic equilib- 
rium arguments. The reported angular rotation Jrequency thresh£)lds_£or vor- 
tex nucleation approximately include 0.1," 0.3,t3 0.4|13 and O.TGJ'EJ'Ej in the 
units of the transverse trap frequency. A number of theoretical investigations 
have been performed to explain these observations.S'iffliS'iiSi'il'fflffl 

The relatively high nucleation thresholds may be qualitatively under- 
stood by considering the energy difference between the vortex and nonyprtex 
states as a function of the position of the vortex in the condensate;Ell in a 
trap rotating at the thermodynamic equilibrium frequency to support an 
axisymmetric vortex line, this self-energy of the vortex features a maximum 
on the condensate surface which prevents the vortex from penetrating the 
system. Thus much higher stirring frequencies are needed in order to remove 
the energy barrier at the condensate surface and hysteretic behavior results. 

In the following, we consider singly quantized, axisymmetric vortex 
states radially confined by a harmonic trapping potential with the frequency 
uj±_ perpendicular to the rotation axis. Although the periodic boundary 
conditions applied in the axial direction makes a direct comparison with a 
specific experimental setup difficult, the relevant physics described by this 
system closest resembles that of oblate/spherical condensate geometries in 
which, for instance, vortex bending is negligible. The axisymmetric con- 
densate states are obtained from the Gross-Pitaevskii equation and the Bo- 
goliubov Eqs. @ are solved for the quasiparticle excitations which serve 
to provide information on the local energetic stabilities of the states. The 
thermodynamic critical angular frequency for vortex stability is determined 
by the condition that the free energies, see Eq. (p!?!), of the vortex and the 
nonvortex states become equal in the rotating frame of reference in which 
the Hamiltonian is time- independent. For technical details of the numerical 
computations performed, see Ref. 
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3.1. Local Energetic Stability 

In nonrotating traps, according to zero-temperature Bogoliubov theory, 
an axisymmetric vortex state is locally unstable in the sense that the vor- 
tex would spiral out of the condensate if the dissipation in the system is 
significant .eHI In other words, the self-energy of the vortex decreases monj> 
tonically as the radial distance of the vortex from the trap center increases .eJ 
Such local energetic instability is also manifest in the spectrum of the vor- 
tex state which contains anomalous mode(s) localized in the vortex core, 
see Fig. |^. These elementary excitations possess negative energies with 
respect to the condensate ground state and thus the condensate could lower 
its energy by transferring particles into the anomalous mode. 

However, when the trap rotation frequency is increased, a plateau in 
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Fig. 1. Computed windows of the Bogoliubov quasiparticle excitation spec- 
tra for a singly quantized vortex state as functions of the angular momentum 
quantum number, /. (a) The spectrum of the nonrotating system exhibits 
one anomalous vortex core localized mode for I = —1 with a negative en- 
ergy with respect to the condensate state. Such an anomalous mode also 
corresponds to the precessing motion of an off-centered vortex line and im- 
plies local energetic instability of the state£^c3 The suggested possible decay 
channel of such a state in the presence of dissipation is the outward spiral- 
ing motion of the vortex leading to vortex annihilation on the condensate 
surface.El (b) The corresponding spectrum computed in the frame rotating 
at an angular frequency close to Qi, see Eq. (18). The anomalous vortex 
core mode {I = —1) is Doppler-shifted to positive energy but the lowest 
surface modes now osculate the condensate energy level. Further increase 
of the rotation would allow for the generation of anomalous surface modes, 
implying an instability of the singly quantized vortex state to a state of an 
increasing number of vortices. 
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the self-energy is formed corresponding to the onset of metastability of the 
axisymmetric vortex state.0 In general, the local energetic stability of the 
condensate state may be identified with the positivity of the quasiparticle 
excitation spectrum of the system. In terms of the Bogoliubov excitation 
spectrum, the local (meta)stability frequency J7m corresponds to the fre- 
quency LOa of the lowest anomalous mode which acquires an energy HQ in 
the frame rotating at the angular frequency $7. 



3.2. Thermodynamic Equilibrium 

As the rotation frequency of the trap exceeds the metastability fre- 
quency Jim; a minimum develops in the self-energy function at the trap 
center. At a sufficiently high frequency, Oc, the free energies of the vortex 
and the nonvortex states evaluated in the rotating frame of reference co- 
incide on the trap axis, indicating the onset of thermodynamic stability of 
an axisymmetric vortex. At this stage, on the condensate boundary there 
still remains the energy barrier preventing the vortex entry which has to be 
overcome. This implies hysteretic behavior when the rotation frequency is 
ramped above the nucleation threshold and back again. 

The (global) thermodynamic stability implies local energetic stability 
of the system. However, in elongated trap geometries the critical frequency 
for the local stability of an axisymmetric vortex line may even exceed the 
corresponding thermodynamic value This is because the true equilib- 

rium configuration in prolate trap geometries corresponds to a bent vortex 
line.iiS The additional rotation required for the local stabilization of an 
axisymmetric vortex line is consumed in the straightening of the otherwise 
bent vortex. 



3.3. Landau Criterion 

In order to explain the high values of the observed vortex nucleation 
thresholds, it has been argued that the correct condition for the removal of 
the energetic barrier on the condensate surface leading to the vortex for- 
mation through the surface instability is given by the Landau criterion for 
the generation of quasiparticle excitations in the systemS Specializing to a 
rotationally invariant systems, the Landau critical angular frequency for the 
creation of the surface modes in the condensate is given by 

0^ = niin{i^}, (18) 
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Fig. 2. Energies per angular momentum for surface modes as functions 
of the angular- momentum quantum number / calculated within the zero- 
temperature Bogoliubov approximation. The solid line is the hydrody- 
namic prediction yielding a fair approximation for the lowest multipole 
modes For higher values of angular momenta, the deviation becomes sig- 
nificant as the hydrodynamic value for loi/1 diverges whereas the Bogoliubov 
spectrum exhibits a minimum corresponding to the Landau critical angular 
frequency, O^. The quadrupole resonance is expected to occur at O « 0.7a;_L 
in agreement with the experimental observations in which the condensate 
has been stirred by a quadrupole drive. 



where hwi is the energy of a quasiparticle carrying the angular momentum hi. 
These surface modes are quasiparticle excitations corresponding to revolv- 
ing density perturbations localized around the condensate surface. Again, in 
terms of the excitation spectrum, the Landau critical frequency corresponds 
to the rotation frequency at which there emerges anomalous negative-energy 
modes localized on the condensate surface in the elementary excitation spec- 
trum of the system, see Fig. ^3. 

However, in a number of recent experiments, the emergence of the vor- 
tices has been seen only for yet higher rotation velocities of the trap. Ac- 
cording to the current understanding, those values are best explained in 
terms of the excitation of the quadrupolar surface modes. In Ref. 11, pat- 
terns of laser beams with different multipolarities were used and the vortex 
nucleation resonances were observed for corresponding frequencies for the 
excitation of multipolar surface modes. Thus it is understood that if the 
rotating anisotropy resonantly excites certain surface modes, the vortex nu- 




Fig. 3. Radial density distributions of the condensate and the nonconden- 
sate components for T ~ Tc/2. For the sake of clarity, the noncondensate 
density is scaled up by a factor of 5. The condensate density vanishes on 
the vortex axis where the noncondensate accumulates forming an additional 
self-stabilizing potential for the vortex state. On the boundary of the con- 
densate, the noncondensate bulges due to the mutual repulsion between the 
two components. 



cleation threshold is determined by the frequency of those excitations rather 
than by the minimum oi ivi/l. 

In the hydrodynamic limit, the dispersion relation for the collective 
quasiparticle excitations of a trapped Bose-Einstein condensate was obtained 
in Ref. For surface modes it reduces to uJi = \flu}^ which is expected to 
be accurate for the lowest angular-momentum states, see Fig. ||. Especially, 
for the quadrupole mode it yields ujijl ~ 0.7li;_l in reasonable agreement 



with the observations reported in Refs. 14,24,25. Although the spectrum. 



and especially depends on the shape of the trapping potential, the value 
for the quadrupoleLxesonance / = 2 is only slightly altered under a change of 



the trap geometry 



28 



Notice also that the thermodynamic critical frequency J7c may be inter- 
preted as the Landau criterion f^c = AF/L for the generation of a vortex 
excitation. Here IS.F = F(vortex) — F(nonvortex) is the 'activation energy' 
needed for the creation of the vortex and L = Nh is its angular momen- 
tum. However, the difference between Vt^ and Vti is that the latter is the 
frequency needed for generating a single excitation whereas the former cor- 
responds to exciting all the particles into the unit angular-momentum state 
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simultaneously. 



3.4. Critical Frequencies at Finite Temperatures 

A characteristic feature for all the self-consistent, first-order finite-tem- 



perature theories considered for vortex states in Refs. |46| and ^ is that none 
of them contain anomalous vortex core modes in distinction to the predic- 
tion of the zero-temperature Bogoliubov approximation. This is because of 
the core-filling effect of the noncondensate component lifting the anomalous 
modes to positive energies, see Fig. |3|. Moreover, it suggests that the ax- 
isymmetric vortex state could be metastable even in the nonrotating trap. 
However, these approximations neglect the noncondensate dynamics causing 
shifts in the energies of the quasiparticle modes. For a further discussion on 
thermal effects on the local stability of vortex states see, for example, Ref. ^ 
and the references therein. 
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Fig. 4. Thermodynamic (•) and Landau (o) critical angular frequencies, 
and the quadrupole resonance (■) for vortex formation as functions of tem- 
perature. All of these frequencies increase monotonically for increasing 
temperature but the temperature dependence is significant only at higher 
temperatures for which the validity of the Popov approximation not treat- 
ing the dynamics of the noncondensate is questionable. The numbers 
above the circles denote the corresponding multipolarity I of the destabiz 
lizing surface mode. The critical temperature for the system is given byt3 
Tc(n) = T^(l — ri^/w^)!/^ where ~ 45 nK for the parameter values used. 
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In Fig. ^ we have plotted the thermodynamic critical frequency (•), the 
critical frequency given by the Landau criterion (o) for generation of surface 
excitations, and the quadrupole resonance (■) as functions of temperature. 
The numbers above the squares denote the multipolarity I of the destabilizing 
surface mode. The decrease in the multipolarity as temperature is increased 
is caused by the shrinking of the condensate with increasing temperature. 
All of the critical frequencies remain fairly constant for temperatures below 
0.6 Tc where the Popov approximation neglecting the off-diagonal mean- 
fields m(r) is expected to be in reasonable agreement with the experiments. 
However, even for higher temperatures the variation of these frequencies is 
such that the critical frequency thresholds for the formation of vortices is 
rather determined by external parameters, such as the geometry of the con- 
fining trap and the shape of the rotating perturbation driving tJie condensate 
towards the vortex state, than the temperature of the system.E3 The results 
for the nucleation thresholds at finite temperatures, see Fig. ^, are in fair 
accordance with the experimentally observed values, see Sec. 3, and with 
those predicted by the zero-temperature field theoriesJHa The quadrupole 
resonance occurs at ~ 0.7 uj±, and for the temperatures accessible in the 
experiments with oblate/spherical condensate shapes for which our compu- 
tational results should closest be applicable, the Landau critical frequency 
yields values close to 0.3u;_l — OAuj±. 



4. DISCUSSION 

The formation process of vortex lattices in a stirred Bose-Einstein cort 
densate studied at MIT exhibited only a weak dependence on temperature.Ea 
Contrary to this, the experiment performed at ENS reported a qualitative 
temperature dependence in the lattice ordering.! In the light of these ex- 
periments and the recent pronounced interest in the vortex formation phe- 
nomena in superfluids, it is instructive to investigate theoretically the finite- 
temperature effects on the vortex nucleation in weakly interacting, gaseous 
Bose-Einstein condensates. 

We have computed the thermodynamic critical frequency Jlc for the 
stability of an axisymmetric vortex line as well as the resonance frequencies 
of the surface modes for a range of temperatures providing, for instance, the 
Landau critical frequency 0,i for vortex nucleation. Only a weak temperature 
dependence in the critical rotation frequencies were found in the temperature 
range within which the Popov approximation applied is expected to yield 
reliable predictions. Thus the conclusion is that the high observed values for 
the vortex nucleation frequency are not solely related to the noncondensate 
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gas. 

It is understood that the lower critical trap rotation frequency Qci for 
the formation of vortices in gaseous Bose-Einstein condensates is rather given 
by the Landau criterion 0,£ for the generation of surface excitations than by 
the thermodynamic value 0,c for the stability of an axisymmetric vortex. If 
the rotating drive of the condensate possesses a well-defined multipolarity /, 
it is the corresponding surface mode which determines the threshold nucle- 
ation frequency, instead of Qg. Moreover, the low stirring frequency reported 
in Ref. 11 is probably explained by the fact that the actual condensate flow 
velocity around the small stirrer beam employed is likely to exceed the rota- 
tional velocity of the beam, hence yielding a much higher effective rotation 
frequency.Ej 

In contrast to the helium superfluids, the limit of an upper critical rota- 
tion, flc2, should be within experimental reach in the gaseous Bose-Einstein 
condensates. Recently, such condensates in rapidly rotating traps have at- 
tracted much interest .tl'ililS When the rotation frequency Q of the trap 
reaches the harmonic trappine frequency uj± , the confinement of the atoms is 
lost and they leave the tra p Close to this limit, the system may enter the 
quantum-Hall-like regimeE3E3E3 and the vortex cores are predicted to melt 
into a single giant multiply quantized vortex. If the atoms could be 

held together above uj±, for instance, by additional stronger than harmonic 
potentials, the physics in the vicinity of the upper critical frequency could 
be experimentally explored. Also the possibility of a Kosterlitz-Thouless 
transition at such extreme rotational conditions has been suggested.E3 

In addition to the giant vortices in the rapidly rotating traps, a gen- 
eration of multiply quantized vorticity in dilute Bose-Einstein condensates 
have been suggested by applying additional pinning potentials for stabiliz- 
ing the multiquantum vortex states £iJ Moreover, vortices with circulation in 
multiplets of two have recently been cr eatedlli by continuously changing the 
direction of the axial magnetic field, as first suggested theoreticallyHIS 
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